Mon. Not. R. Astron. Soc. 000, dlSl (2005) Printed 2 February 2008 (MN WF^ style file v2.2) 



The open cluster initial-final mass relationship and the 
high-mass tail of the white dwarf distribution 

Lilia Ferrario ^, Dayal Wickramasinghe^, James Liebert^ and Kurtis A. Williams^ 

^ The Australian National University, ACT 0200, Australia 

'^Steward Observatory, University of Arizona, Tucson, AZ 85726, USA 



ABSTRACT 

Recent studies of white dwarfs in open clusters have provided new constraints on the 
initial - final mass relationship (IFMR) for main sequence stars with masses in the 
range 2.5 — 6.5Mq. We re-evaluate the ensemble of data that determines the IFMR and 
argue that the IFMR can be characterised by a mean initial-final mass relationship 
about which there is an intrinsic scatter. We investigate the consequences of the IFMR 
for the observed mass distribution of field white dwarfs using population synthesis 
calculations. We show that while a linear IFMR predicts a mass distribution that is in 
reasonable agreement with the recent results from the PG survey, the data are better 
fitted by an IFMR with some curvature. Our calculations indicate that a significant 
(~ 28) percentage of white dwarfs originating from single star evolution have masses in 
excess of ~ 0.8Mq obviating the necessity for postulating the existence of a dominant 
population of high-mass white dwarfs that arise from binary star mergers. 

Key words: Stars: white dwarfs - stars: evolution; open clusters and associa- 
tions: individual: Pleiades, NGC2516, NGC2168 (M35), Praesepe, Hyades, NGC3532, 
NGC2099 (M37), Sirius. 



1 INTRODUCTION 

A knowledge of the Initial - Final mass relationship (IFMR) 
of stars is of fundamental importance in many astrophysi- 
cal contexts and plays a key role in our understanding of 
chemical enrichment and the efficiency of star formation in 
galaxies. The low and intermediate mass end of this rela- 
tionship can be probed by studying the properties of white 
dwarfs in open clusters. 

The IFMR can also be investigated by looking at the 
mass distribution of field white dwarfs. Here, we are dealing 
with white dwarfs which would have had an origin in stel- 
lar clusters with a variety of parameters during the history 
of the galactic disc. Earlier studies (e.g. Koester & Weide- 
mann 1980) showed that the predominant characteristic of 
the mass distribution, a sharp peak at Mj = 0.57±.01, could 
be explained by the steep initial mass function and the finite 
age of the galactic disc. This peak corresponded very closely 
to the current turn-off mass of the stellar population in the 
galactic disc. The nature of the mass distribution above and 
below this peak clearly contained important additional in- 
formation on the progenitors of the white dwarfs, but this 
could not be fully exploited at that time because of lack of 
data. 

Recently, detailed studies have been carried out on large 
samples of white dwarfs in surveys such as the PG and the 



SDSS. With these new data, the relative importance of mas- 
sive (A/ > 0.8Mq) white dwarfs only became apparent when 
attempts were made to correct observations based on mag- 
nitude limited surveys for the mass dependence of the in- 
trinsic luminosity of the white dwarfs. The first such results 
for the PG survey were presented for the magnetic white 
dwarfs by Liebert Bergeron & Holberg (2003). They showed 
that when corrections were made to allow for the fact that 
the more massive magnetic white dwarfs had smaller radii, 
and hence were less luminous at a given effective tempera- 
ture, the space density of high-field magnetic white dwarfs 
(HFMWDs) more than doubled from the previously esti- 
mated value of ~ 4 percent to ~ 10 percent. Wickramas- 
inghe & Ferrario (2005, hereafter WF05) used population 
synthesis calculations to show that these observations could 
be explained if a significant fraction of the HFMWDs had 
massive progenitors. The mass distribution of the entire DA 
white dwarf population (magnetic and non-magnetic) in the 
PG sample has only recently been established by Liebert, 
Bergeron & Holberg (2005a). These data, when corrected to 
an equivalent volume limited sample, showed the presence of 
a significant high-mass tail with ~ 30 per cent of DA white 
dwarfs having masses in excess of O.8A/0. The details of this 
mass distribution must hold important clues on the nature 
of the mean IFMR appropriate to stars in the galactic disc. 

In this paper we re-evaluate the open cluster data that 
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are currently being used to define an IFMR and characterise 
this relationship. We show that the nature of this relation- 
ship has observable consequences for the mass distribution 
of field white dwarfs and show how the latter can be used 
to place constraints on the IFMR. 



2 ANALYSIS OF CLUSTER DATA 

We use data for the clusters listed below. Only white dwarfs 
which are probable cluster members are included from the 
published data sets. The masses of the progenitor stars are 
calculated using the nuclear lifetimes from the tables of Gi- 
rardi et al. (2002) for Z = 0.008, 0.019, 0.03. The white dwarf 
mass and cooling age are obtained from the Fontaine, Bras- 
sard, & Bergeron (2001) evolutionary models which incorpo- 
rate improvements in the equation of state and conductive 
opacities. 

Hyades and Praesepe — The data points are from 
Claver et al. (2001) supplemented with recent results on two 
new stars in Praesepe (Dobbie et al. 2005). Ferryman et al. 
(1998) estimate an age 625jl5Q Myr with the errors given 
by von Hippcl (2005) for the Hyades. The isoclirone fit in 
Claver et al. (2001) led these authors to adopt the same age 
for Praesepe. Both have been linked to the same Hyades 
Supcrclustcr and have a metallicity of [Fe/H] — 0.11 — 0.15 
(King & Hiltgcn 1996; Cayrel, Cayrel de Strobel & Campbell 
1985; Bocsgaard & Friel 1990). We have used Z = 0.019. 

Pleiades — The age of the Pleiades is I25I25 from 
isochrones calculated with moderate convective core over- 
shooting (Mazzci & Pigatto 1989; Meynet, Mermilliod & 
Maeder 1993; Ventura ct al. 1998). This is in agreement with 
the age determined from the luminosity at which lithium is 
detected in the substellar objects (Stauffer, Schulz, & Kirk- 
patrick 1998). Since the metallicity of the Pleiades does not 
differ more than ±0.1 dcx from solar (e.g. King & Sodcrblow 
2000), the Girardi et al. (2002) Z = 0.019 table was used. 
We have adopted the T^ff obtained from the spectral fit, but 
have modified the gravity to give the gravitational redshift 
mass of I.O2M0. 

NGC2168 (M35) - The cluster data are from Williams 
et al. (2004) supplemented by new data on LAWDS 2 and 
LAWDS 27. We adopt an age 150 ± 60 Myr given by von 
Hippel (2005) and used by Williams et al. (2004). A separate 
study by Steinhauer (2003) gives an age of 160 ± 20 Myr. 
Sung & BcsscU (1999) obtain a value of [Fc/H] = -0.3 while 
Twarog et al. (1997) give [Fe/H] = -0.16. In our studies, we 
use Z = 0.008. 

NGC2516 — The cluster data are from Koester & 
Reimers (1996). Meynet et al. 1993 estimate an age 158^20 
Myr. Sung et al. (2002) give an abundance [Fe/H]= -0.10± 
0.04. We have used Z = 0.019. 

NGC3532 — The cluster data are from Koester & 
Reimers (1993). The \ogg and T^jj used by us are the mean 
values of the observations in Koester & Reimers (1996). This 
sparse cluster has an estimated age of 30015^50 M^yi'- Twarog 
et al. (1997) give a metallicity [Fe/H] = -0.022. Thus, we 
have used the Z = 0.019 table. 

NGC2099 (M37) Wc have used the results of Kali- 
rai et al. (2005). This set of data are of somewhat lower 
quality because these stars are much fainter than all other 
stars considered here (F ~ 22 — 23.6) and as a result, the 



coverage of the high Balmer lines which are most sensitive 
to gravity (and stellar mass) was generally of poorer qual- 
ity. In addition, the stars WD 1, 15, 17, 24 are excluded 
because they are at the wrong distance to be members of 
this cluster. We note also that WD 15 and 17 have cooling 
ages that are much larger than any plausible cluster age. 
WD6 and WD21 are also excluded because the gravity was 
determined assuming cluster membership, rather than from 
the spectrum. We have adopted the age of 520tso ^^'^ solar 
metaUicity (Twarog ct al. 1997). 

Sirius — We include new data by Liebert et al. (2005b) 
who used the well determined luminosity and radius of Sirius 
A and stellar evolution models with a solar metallicity to 
accurately constrain the age for this system to 238!t}3. 

Error bars for the initial and final masses were deter- 
mined using a Monte Carlo simulation. 1000 white dwarfs 
were created with Gaussian distributions of Teff and logg 
with 1(7 equivalent to the published errors on these quan- 
tities. This makes the implicit assumption that the derived 
Teff and \ogg are independent. The Mi and Mj were then 
calculated for each simulated sample, and the 68 percent 
confidence limits. These results are presented together with 
the white dwarf cooling ages r in Table 1. Two sets of results 
are given for the errors. The first assumes the best estimate 
for the age of a cluster so that the errors are those from the 
white dwarf spectral observations. The second gives the sys- 
tematic errors that result purely from the uncertainties on 
the cluster ages. The total error bars are obtained by adding 
these errors in quadrature as shown in Figure 1. 



3 INITIAL TO FINAL MASS RELATIONSHIP 

The data in Figure 1 span the initial mass range ~ 2.5 — 
6.5Mq and thuse do not provide information at the low and 
high mass ends of single stars evolving into white dwarfs. 
We note that the IFMR of most clusters appears to have 
an intrinsic spread, that is, there is a range of initial masses 
which gives rise to a given final mass. If this spread is real 
(that is, if we can dismiss the possibility that the spread is 
not caused, for example, by non membership of the cluster, 
effects introduced by binaries, or by the underestimation of 
errors, as could be the case for NGC2099) then the IFMR 
should be characterised by a mean relationship and a spread 
about this mean. This spread could in principle also be mod- 
elled in future when enough data become available. 

Given the various uncertainties involved, we cannot jus- 
tify the use of any but a linear relationship to model the 
cluster data. When allowance is made for the uncertainties 
in the cluster ages, the error bars increase significantly par- 
ticularly at the high mass end, and the weight is shifted 
from the high-mass end to the low-mass end of this dia- 
gram. Given the scatter in this part of the diagram, a least 
squares linear fit to the data leads to a meaningless result. 
We have circumvented this problem by introducing an an- 
chor point at (Mi = 1.1 ± O.IMq.M/ = 0.55 ± O.OIMq) 
(c.f. Weidemanri 2000) which we justify as a requirement to 
reproduce the well established peak in the field white dwarf 
mass distribution at M/ = 0.57Mo (see next section). With 
this anchor point, we obtain the best fit linear relationship 
to be 

Mf = (0.10038 ± 0.00518) Mi + 0.43443 ± 0.01467 
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Table 1. Cluster data. 



Cluster 


White Dwarf 


Teff ± dT^ff 


logg ± dlogg 


Mf ± dMf 


log T ± d log r 


Mi 


dMi 


dMi 






















(random) 


(systematic) 


Pleiades 


LB 1497 


31660 ± 350 


8.64 




l.Uzo ± U. 


.026 


7.781 ± 0.062 


6. 


,542 


+ 0.458 
— 0.346 


+ 1, 

-0 




Haydes 


0352+098 


16629 ± 350 


8.16 + 


.05 


u. i iy ± u. 


030 


8.270 ± 0.046 


.3. 


094 


+ 0.052 
— 0.045 


+ 0, 
-0 


'l34 
,114 


& 


0406+169 


15180 ± 350 


8.30 ± 


.05 


U.oUb ± U. 


.031 


8.488 ± 0.046 


3. 


165 


+0.147 
— 0.109 


+ 0, 
-0 


22-:} 
,172 


Praesepe 


0421+162 


19570 ± 350 


8.09 ± 


.05 


U.D80 ± 0. 


.031 


7.970 ± 0.055 


2. 


892 


+0.023 
— 0.020 


+ 0, 

-0 


103 
,090 




0425+168 


24420 ± 350 


8.11 ± 0. 


.05 


U.7U5 ± 0. 


.031 


7.549 ± 0.077 


2. 


.789 


+ 0.010 
— 0.011 


+ 0, 
-0 


()88 
,079 




0431+125 


21340 ± 350 


8.04 ± 0. 


.05 


U.D52 ± 0. 


.032 


7.752 ± 0.068 


2. 


.825 


+ 0.017 
— 0.014 


+ 0, 
-0 


()93 
,083 




0438+108 


27390 ± 350 


8.07 ± 


.05 


U.D84 ± 0. 


.031 


7.203 ± 0.060 


2. 


.757 


+ 0.004 
— 0.003 


+ 0, 
-0 


()84 
,076 




0437+138 


15335 ± 350 


8.26 + 0, 


.05 


U.782 ± 0. 


.033 


8.447 ± 0.047 


3. 


,365 


+ 0.117 
— 0.089 


+ 0, 
-0 


192 
,156 




0836+197 


21949 ± 350 


8.45 + 0, 


.05 


n nnri i r\ 

u.9uy ± u. 


.030 


8.136 ± 0.049 


2. 


,981 


+0.035 
-0.031 


+0, 

-0, 


115 
,101 




0836+201 


16629 ± 350 


8.01+0. 


.05 


U.DzU ± U. 


.031 


8.154 ± 0.050 


2. 


.993 


+0.046 
— 0.028 


+ 0, 

-0 


117 
,102 




0836+199 


14060 ± 630 


8.34 + 0. 


.06 


n Q01 _|_ n 


.uo / 


8.612 ± 0.069 


Q 
O. 


QQ7 


+ 0.652 
— 0.350 


+ 0, 
-0 


,297 




0837+199 


17098 ± 350 


8.32 + 0, 


.05 


n Qi n 1 n 
U.oiy ±: U. 


.Uoz 


8.351 + 0.048 


Q 
O. 




+ 0.081 


+0, 

-0, 


153 
,129 




0837+218 


17845 ± 560 


8.48 + 0, 


.075 


u.yio ± u. 


.045 


8.422 ± 0.073 


3. 


.315 


+ 0:i68 
—0. 124 


+0, 

-0, 


178 
,148 




0837+185 


14170 ± 1485 


8.46 + 0, 


.15 


u.yuu ± U. 


.088 


8.683 ± 0.171 


4. 


.695 


+9.990 
— 1.131 


+0, 

-0 


901 

,523 




0840+200 


14178 ± 350 


8.23 ± 0. 


.05 


0.761 di 0. 


.033 


8.522 ± 0.045 


3. 


.572 


+0.186 
—0.129 


+ 0, 
-0 


255 
,196 


NGC2516 


2516-1 


28170 ± 310 


8.48 ± 0. 


.17 


n QQ1 -U n 
u.yo-L ^ u. 


HQS 


7.760 ± 0.230 


u. 


411 


+ 1.017 
— 0.500 


+ 0, 

-0 


524 
,386 




2516-2 


34200 ± 610 


8.60 + 0, 


.11 




.uuo 


7.621 + 0.179 


c; 
o. 


.uyu 


+ 0.406 


+ 0, 
-0 


,311 




2516-3 


26870 ± 330 


8.55 + 0, 


.07 


n c\ci.c\ _L r\ 


.039 


7.922 + 0.079 


6. 


,141 


+ o!692 


+ 0, 
-0 


924 
,587 




2516-5 


30760 ± 420 


8.70 + 0, 


.12 


1.054 i: 0. 


.063 


7.883 ± 0.136 


5. 


.902 


+ a986 
—0.524 


+0, 

-0 


786 

,514 


NGC3532 


3532-8 


23367 ± 1065 


7.713 + 0.148 


U.4oi ± U. 


.065 


7.324 ± 0.098 


3. 


.634 


+ 0.028 
— 0.014 


+ 1 , 
-0 


259 

,514 




3532-9 


29800 ± 616 


7.827 ± 0.229 


n cci _|_ n 


111 


6.960 + 0.107 


Q 
O. 


f o 


+ 0.008 
— 0.002 


+ 1, 
-0 


142 

,488 




3532-10 


19267 ± 974 


8.143 + 0.266 


U.712 ± 0. 


.158 


8.043 ± 0.244 


4. 


.200 


+ 0.937 


+ 3, 
-0 


966 

,815 


NGC2099 


2099- WD2 


19900 ± 900 


8.11 + 0, 


.16 


0.690 ± 0. 


.098 


7.961 ± 0.165 


3. 


,120 


+o!io5 

— 0.071 


+ 0, 

-0, 


,162 


(M37) 


2099- WD3 


18300 ± 900 


8.23 + 0, 


.21 


0. m4 ± 0. 


,128 


8.190 ± 0.190 


3. 


.300 


+0.282 
—0.162 


+0, 
-0, 


258 
,197 




2099- WD4 


16900 ± 1100 


8.40 ± 0. 


.26 


0.870 ± 0. 


.153 


8.428 + 0.217 


3. 


.773 


+2.146 
— 0.435 


+0, 

-0 


495 

,324 
305 
,219 




2099- WD5 


18300 ± 1000 


8.33 ± 0. 


.22 




1 '^1 

. AO A 


8.271 ± 0.188 


Q 
O. 


.^uu 


+ 0.453 
— 0.200 


+0, 

-0 




2099- WD7 


17800 ± 1400 


8.42 + 0, 


.32 


U.883 ± 0. 


.196 


8.378 ± 0.289 


3. 


.623 


+ 2.733 


+ 0, 
-0 


405 
,278 




2099-WD9 


15300 ± 400 


8.00 + 0, 


.08 


U.D13 ± 0. 


.053 


8.268 ± 0.067 


3. 


.401 


+ 0!l22 
— 0.087 


+ 0, 
-0 


302 
,218 




2099- WDIO 


19300 ± 400 


8.20 + 0, 


.07 


0.748 + 0. 


,043 


8.090 ± 0.068 


3. 


,205 


+0.065 
—0.044 


+0, 

-0, 


225 
,178 




2099- WD 11 


23000 ± 600 


8.54 ± 0. 


.10 


0.961 ± 0. 


.056 


8.149 + 0.095 


3. 


.257 


+0.104 
-0.079 


+0, 

-0 


242 
,188 




2099- WD12 


13300 ± 1000 


7.91+0. 


.12 


0.550 ± 0. 


.071 


8.394 + 0.126 


3. 


,665 


+ 0.575 
-0.252 


+ 0, 
-0 


429 
,290 




2099- WD13 


18200 ± 400 


8.27 + 0, 


.08 


0.789 ± 0. 


,052 


8.229 + 0.074 


3. 


,348 


+ 0.110 
-0.085 


+ 0, 
-0 


278 
,207 




2099- WD14 


11400 + 200 


7.73 + 0, 


.16 


0.450 + 0. 


,079 


8.487 + 0.085 


4. 


,011 


+ 0.623 
-0.326 


+ 0, 

-0, 


674 
,400 




2099- WD16 


13100 ± 500 


8.34 + 0, 


.10 


0.825 + 0. 


,062 


8.689 + 0.082 


9. 


,136 


+9.990 
-3.864 


+9, 

-3, 


,990 
,737 


NGC2168 


NGC2168 LAWDS 1 


32400 ± 512 


8.40 ± 0. 


.125 


0.894 ± 0. 


.074 


7.371 ± 0.220 


4. 


.839 


+0.303 
-0.163 


+1, 

-0 


642 
,750 


(M35) 


NGC2168 LAWDS 2 


32700 ± 603 


8.34 + 0. 


.080 


0.859 ± 0. 


,048 


7.232 ± 0.148 


4. 


,735 


+ 0.123 
-0.077 


+1, 

-0 


464 
,704 




NGC2168 LAWDS 5 


52600 ± 1160 


8.24 + 0, 


.095 


0.828 + 0. 


,053 


6.223 + 0.047 


4. 


,508 


+ 0.002 
-0.003 


+1, 

-0 


164 
,602 




NGC2168 LAWDS 6 


55200 ± 897 


8.28 + 0, 


.065 


0.852 + 0. 


,038 


6.130 + 0.048 


4. 


.504 


+ 0.002 
-0.002 


+1, 

-0, 


159 
,601 




NGC2168 LAWDS 15 


29900 ± 318 


8.48 + 0, 


.060 


0.934 + 0. 


,034 


7.652 ± 0.093 


5. 


,239 


+0.238 
-0.192 


+2, 

-0, 


,618 
,932 




NGC2168 LAWDS 27 


30500 ± 397 


8.52 ± 0. 


.061 


0.957 ± 0. 


.035 


7.676 ± 0.092 


5. 


.298 


+0.260 
-0.217 


+2, 
-0 


793 
,963 


Sirius 


Sirius B 


25000 ± 200 


8.60 ± 0. 


.04 


1.000 ± 0. 


.020 


8.091 ± 0.020 


5. 


.056 


+0.262 
-0.171 


+ 0, 
-0, 


273 
,213 



We note that if we adopt the best estimate for the 
cluster ages, a best fit Unear relationship can be obtained 
from the cluster data alone without an anchor point at the 
low-mass end. The linear relationship that we obtain with 
this assumption is indistinguishable from the above equation 
within the given errors. 



4 POPULATION SYNTHESIS 

We show in Figure 2 the expected mass distribution of DA 
white dwarfs deduced from the PG survey corrected to an 
equivalent volume limited sample (Liebert et al. 2005a). 

Caxbon-Oxygen white dwarfs resulting from single star 
evolution are not expected to have a mass less than O.SM©. 



The white dwarfs with lower masses in the PG survey must 
have interiors composed of helium and result from binary 
evolution during which the envelopes have been stripped off 
before helium ignition. We do not consider such stars in our 
population synthesis calculations. 

We have carried out calculations of the expected mass 
distribution of field white dwarfs for different assumptions 
on the IFMR. We follow closely the method described by 
WF05 and we refer to that paper for details. 

In essence, we calculate the number of white dwarfs per 
unit volume per unit final mass interval with luminosities in 
the range Li to L2 by the integral 



dNwDiMf,Li,L2) 
dMf 



Jtc H{Tdisc — 



— Tv0(MO^dt 
t. Mi) ' 'dMf 
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Figure 1. The ensemble of data for the seven open clusters 
defining the initial to final mass relation. The error bars cor- 
respond to the best estimates of the cluster ages. Red: Pleiades; 
Cyan: NGC2516; Magenta: NGC2168 (M35); Blue: Praesepe and 
Hyades; Green: NGC3532; Yellow: NGC2099 (M37); Black: Sir- 
ius. The solid line is the best linear fit to the data and the dashed 
line is the IFMR with curvature (see text). 

with: 

Tc = Tdisc — Tevol{Mi) — Tu,d{M f (Mi) , L2) 



T 0.2 



1 1 1 1 




1 


1 1 1 1 1 1 1 1 1 1 1 1 
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Figure 2. Weighted observed mass distribution of field white 
dwarfs in the Palomar-Green survey (Liebert et al. 2005a, shaded 
histogram) for stars hotter than 13,000 K with the best linear 
IFMR fit model (top panel) and the IFMR with curvature (bot- 
tom panel). White dwarfs with masses less than 0.5Mq are be- 
lieved to be helium white dwarfs resulting from close binary evo- 
lution and are not modelled in this paper. 



Td = Tdisc - Tevol{Mi) - T^d{Mf{Mi), Ll) 

and where 4>{Mi) is the initial mass function, tl;{t) is the star 
formation rate per unit area of the galactic disc, Tdisc is the 
age of the disc, Tevoi is the time of evolution of a Main Se- 
quence star of mass Mi to the beginning of the white dwarf 
phase (we assume that all white dwarfs are born at luminos- 
ity lO^Z/0), t^d is the white dwarf cooling time to a luminos- 
ity L and H {Tdisc — t, Mi) is a possibly time-dependent scale 
factor which allows for the inflation of the galactic disc. The 
parameters used here are as in WF05. 

At the basis of our calculations is the assumed form 
for the IFMR. We consider only those stars with Li > 
2.5 Z/0, which corresponds to the luminosity of a white 
dwarf of M = 0.8Mq at Teg = 13,000 (Wood 1995). We 
consider this to be appropriate in an average sense for com- 
parison of our results with the data histogram obtained from 
the PG survey. We show in Figure 2 the results that we ob- 
tain when we use the linear IFMR overlapped to the data. 
The theoretical predictions are in overall agreement with 
the observations for M > O.5M0 although the shape is not 
reproduced in detail. The model predicts ~ 29 percent of 
white dwarfs of all types to have masses greater than O.SMq. 
This value is similar to the ~ 30 percent in the PG survey 
(Liebert et al 2005a). The inclusion of non-DA spectral types 
could also have a small (~ 10 percent) effect on the relative 
numbers and percentages. 

The observed white dwarf mass distribution appears to 
show a deflcit of stars relative to the model in the mass range 
0.65 — O.SMq and an excess in the mass range 1.0 — IAMq. 
Rapid changes in the slope of the IFMR are reflected in the 
white dwarf mass distribution so that these deviations can 
be explained if we move away from a linear IFMR. A well 



determined white dwarf mass distribution can in principle be 
inverted to establish the appropriate IFMR. We demonstrate 
that the introduction of curvature into the IFMR in the 
appropriate mass ranges does indeed result in better overall 
agreement with the observed white dwarf mass distribution. 
The relationship that we have used is given by 

M/ = - 0.00012336Mf + 0.003160Aff - 0.02960M* 

+ 0.12350Mf - 0.21550M^ + 0.19022A/i + 0.46575 

With this relationship, there are changes in slope in the 
white dwarf mass distribution at O.65M0 and at O.95M0 as 
indicated by the PG survey. The cluster data may provide 
independent evidence for the presence of such a curvature 
when a larger and better refined dataset becomes available 
through future cluster observations. 

The interpretation of the mass distribution of field white 
dwarfs is complicated by the fact that there will be a com- 
ponent that has resulted from close binary evolution. We 
expect in particular that white dwarf mergers will play a 
dominant role at the very high mass end (M^ IMq) of this 
distribution. 



5 DISCUSSION AND CONCLUSIONS 

The initial-final mass distribution can be derived empirically 
by studying white dwarfs in open clusters, or by modelling 
the observed mass distribution of white dwarfs. Uncertain- 
ties in cluster ages, metallicity, and in atmospheric model 
fitting to spectral data result in larger error bars being as- 
sociated with many of the cluster data points. In this paper, 
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we have brought together all presently available data on the 
IFMR and made the simplest hypothesis that the the IFMR 
can be modelled by a mean relationship about which there is 
an intrinsic scatter. We have taken the point of view that a 
linear fit is all that can be warranted by the presently avail- 
able cluster data and that the data set is not sufficiently 
extensive to warrant the modelling of the scatter. 

The IFMR can be used to predict the white dwarf mass 
distribution. Here, the analysis depends on the input to the 
population synthesis calculations, such as the well estab- 
lished initial mass function, and the age of the galactic disc. 
The peak in the observed mass distribution is determined by 
the lowest mass main sequence stars that evolve into white 
dwarfs, while the high-mass tail is determined by the pro- 
file of the IFMR for high-mass progenitors. We have shown 
while the best fit linear relationship to the ensemble of clus- 
ter data with a suitably chosen anchor point at the low mass 
end predicts a white dwarf mass distribution that is in gen- 
eral agreement with observations, a relationship with curva- 
ture improves the quality of the fit to the white dwarf mass 
distribution. For both the linear and curved IFMRs, some 
~ 28 percent of the stars in the white dwarf mass distribu- 
tion has masses above ~ O.SMq. It is therefore not necessary 
to postulate the presence of a dominant component in the 
high-mass tail that arises from mergers (e.g. Liebert et al. 
2005a), even though we may expect more mergers in high 
galactic latitude surveys such as the PG, compared to, for 
example, the EUV survey which roughly tracks the B-type 
stars of Gould's belt. 

An interesting corollary to our present study relates to 
the properties of the HFMWDs. It has been know for quite 
some time that the HFMWDs had masses that were, on 
the average, higher than for the non-magnetic white dwarfs. 
Early studies suggested that ~ 4 percent of the white dwarfs 
were magnetic, and such a percentage could be explained 
by assuming that their progenitors were the known mag- 
netic Ap and Bp stars which comprised ~ 10 percent of all 
A and B stars. However, when appropriate corrections were 
made to allow for the fact that most surveys were magnitude 
limited, the predicted proportion of magnetic white dwarfs 
increased by a factor of ~ 2 — 3 (Kawka 2004, Liebert et 
al. 2003). WF05 argued that if one dismissed the possibil- 
ity of a IFMR that is affected by magnetic fields, both the 
observed mass distribution, and the observed space density 
of the HFMWDs could be reproduced if in addition to the 
currently recognised group of magnetic Ap and Bp stars, 40 
percent of stars with masses in the range 4.5 — WMq also 
evolved into HFMWDs. 

This conclusion is now seen to be consistent with the 
results of the present study where 40 — 50 percent of the 
~ 28 percent of stars in the high-mass tail of the white 
dwarf mciss distribution provides the observed ~ 10 percent 
of the HFMWDs. 
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